Constraints from primordial black hole formation at the end of inflation 
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Primordial black hole (PBH) abundance limits constrain the primordial power spectrum, and 
hence models of inflation, on scales far smaller than those probed by cosmological observations. 
Single field inflation models which are compatible with all cosmological data can have large enough 
perturbations on small scales to overproduce PBHs, and hence be excluded. The standard formulae 
for the amplitude of perturbations do not hold for modes that exit the horizon close to the end of 
inflation however. We use a modified flow analysis to identify models of inflation where the amplitude 
of perturbations on small scales is large. For these models we then carry out a numerical evolution 
of the perturbations and use the PBH constraints on the power spectrum to eliminate models which 
overproduce PBHs. Significant PBH formation can occur in models in which inflation can continue 
indefinitely and is ended via a secondary mechanism. We demonstrate that PBHs constrain these 
types of inflation models and show that a numerical evaluation of the power spectrum decreases the 
number of otherwise viable models of inflation. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Primordial black holes (PBHs) can form in the early 
Universe via the collapse of large density perturba- 
tions 1,2]. There are tight constraints on the abundance 
of PBHs formed from their present day gravitational ef- 
fects and the consequences of their evaporation. These 
limits can be used to constrain the power spectrum of 
the primordial density, or curvature, perturbations. The 
PBH constraints on the curvature power spectrum are 
fairly weak, being many orders of magnitude larger than 
the measurements on cosmological scales. They do, how- 
ever, apply over a very wide range of scales and therefore 
provide a useful constraint on models of inflation. Peiris 
and Easther Q have shown that there are single field in- 
flation models, which are compatible with all cosmolog- 
ical observations, for which the perturbation amplitude 
on small scales is large enough to produce a significant 
density of PBHs. 

For scales which exit the horizon close to the end of in- 
flation the standard (Stewart-Lyth [3]) formulae for the 
amplitude of perturbations do not hold. Leach and Lid- 
die Q carried out a numerical calculation of the evolution 
of perturbations for a quadratic inflationary potential. 
They found that the perturbations on scales which exit 
the horizon close to the end of inflation were roughly an 
order of magnitude larger than predicted by the Stewart- 
Lyth formula (see also Ref . |6[ ) • Therefore to fully exploit 
the power of PBH constraints on inflation models, a nu- 
merical calculation of the amplitude of perturbations on 
small scales is required. It has recently been shown 0, H| 
that PBHs can also form on scales which never leave the 
horizon. We do not consider this possibility here. 



In this paper we use a modified flow analysis to identify 
inflation models where the perturbations at the end of 
inflation may be large enough for primordial black holes 
to be overproduced. For these models we carry out a 
numerical evolution of the primordial perturbations and 
use the PBH constraints on the power spectrum to elim- 
inate models which overproduce PBHs. We describe the 
modified flow analysis in Sec. Ill Al and the evolution of 
perturbations and the calculation of the power spectrum 
in Sec. Ill Bl We apply the primordial black hole abun- 
dance constraints and present our results in Sec. IHII and 
conclude with discussion in Sec. IIVI 



II. METHOD 
A. Flow equations approach 

We consider the Hubble slow roll- parameters [|| : 
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where m-p\ is the Planck mass and ' denotes differentia- 
tion with respect to the scalar field, 4>. The flow equa- 
tions [l(| [H| encode the variation of the slow- roll param- 
eters in terms of the number of e-foldings from the end 
of inflation, N = In [a(t cn d) / a(t)] and provide a method 
for stochastically generating inflation models: 
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where oh = 2{ 1 \h) — 4e# and 1 Xh = i]h- 

Following Kinney [Tl| we randomly chose 'initial' val- 
ues for the slow- roll parameters and N cos , the number of 
e-foldings between cosmological scales exiting the horizon 
and the end of inflation, in the ranges: 
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[40,60], 
[0,0.8], 
[-0.5,0.5] , 
[-0.05,0.05] , 
[-0.005,0.005] 
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truncating the hierarchy at M = 6. We then evolve the 
flow equations forward in time (diV < 0) from N = N cos 
until either N — or inflation ends with ejj = 1. In the 
former case we calculate the cosmological observables, 
the spectral index, n s , its running, dn s /dlnfc, and the 
scalar to tensor ratio, r, using the initial values of the 
slow- roll parameters [111 ]: 
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where d = 4(ln2 + 7) - 5 « 0.0814514 and 7 « 0.577. 
In the latter case we evolve the flow equations backward 
Nqos e-folds and calculate the cosmological observables at 
this point. In some cases inflation also ends when evolv- 
ing backwards before N cos e-folds are achieved. These 
models are incapable of supporting the required amount 
of inflation and are discarded. 

Our algorithm differs from that originally proposed by 
Kinney [Tljj in how we handle models chosen from the ini- 
tial hierarchy that are destined to inflate forever, e# — » 0, 
but do not reach this limit within N cos e-foldings. In 
the original flow algorithm in this case the cosmological 
observables are calculated at the late-time fixed point 
i.e. the model is forced to evolve to its asymptotic limit. 
In this limit the running of the spectral index is negli- 
gible. Therefore for models which are compatible with 
the WMAP 7 year measurement of the spectral index, 
n s = 0.964 ± 0.012 [12J, the amplitude of the curvature 
perturbations can not be large on an y sc ale and PBHs arc 
never formed in significant numbers [13j . Following Peiris 
and Easther 0] , we do not force these models to evolve to 
their asymptotic limit but instead terminate them once 
N cos e-folds of inflation have occurred. At this point it is 



assumed that another mechanism, for example a second- 
field such as in hybrid inflation [blj , terminates inflation. 
With this treatment some of these models are consistent 
with the WMAP measurements of the spectral index and 
its running, but have perturbations on small scales which 
may be large enough to over-produce PBHs [3, fT5T - li~7| . 



B. Perturbation calculation 

The evolution of inflationary curvature perturbations, 
1Z, is carried out using the Mukhanov variable [llj], u = 
—zTZ, where 
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The Fourier modes, Ufc, evolve according to a Klein- 
Gordon equation with a time-varying effective mass: 
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where r is conformal time, dr = dt/a, and 
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At early times, ti, when a mode k is well within the 
horizon, aH/k 0, the initial condition for Uk{ri), is 
taken to be the Bunch-Davies vacuum state, 



u k {n) 



1 



exp (— ikr\) . 



(13) 



In the super horizon limit, k 2 <C z", eq. (|11[) has a growing 
mode solution Uk oc z, so that the curvature perturbation 
Tlk = \uk/z\ 'freezes out' and becomes constant. The 
power-spectrum of the curvature perturbations can thus 
be calculated as 
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Eq. (fTT]) can be solved exactly for the special case of 
power-law inflation. The commonly used Stewart-Lyth 
formula is found via a slow-roll expansion around this 
exact solution 14]: 
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where C = —2 + In 2 + 7 w —0.729. This expression 
gives the power spectrum in the asymptotic supcrhori- 
zon limit, k/aH — > 0, in terms of the Hubble parame- 
ter and slow-roll parameters evaluated at horizon cross- 
ing [l9| . It is valid provided that the slow-roll approxima- 
tion holds (specifically that the slow-roll parameters are 
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FIG. 1: The power spectrum of the primordial curvature per- 
turbations generated during the final few e-folds for an exam- 
ple inflation model. The black solid line shows the Stewart- 
Lyth equation while the blue dotted line is the result of a 
numerical mode by mode evolution. 



FIG. 2: The cosmological observables, n s , dn s /dlnfc and r, for 
models generated using the modified flow algorithm described 
in the text. All models which sustain the required number of 
e-foldings of inflation are shown in the top two plots and those 
remaining once PBH constraints are applied are shown in the 
bottom plots. 



slowly varying around horizon crossing) and the asymp- 
totic limit is reached before inflation ends [2(| • For modes 
which exit the horizon close to the end of inflation the 
asymptotic limit will not be reached, and the slow-roll 
approximation may also be violated. Leach & Liddle [|[ 
investigated this for a simple quadratic chaotic inflation 
model. They found that for scales that exit the horizon 
very close to the end of inflation the power spectrum is 
roughly an order of magnitude larger than that found us- 
ing the Stewart-Lyth expression. In other words, analytic 
calculations can significantly underestimate the ampli- 
tude of perturbations and hence the abundance of PBHs 
formed. Therefore, a numerical calculation of the pertur- 
bation evolution is required to accurately compute the 
primordial power spectrum on the very smallest scales. 

We use a modified version of the Inflation v2 module 
(written by Lesgourgues & Valkenburg) [2l[ to carry out 
an accurate numerical calculation of the evolution of per- 
turbations. Fig. [T]shows the power spectrum of curvature 
perturbations of an example inflation model generated 
using the modified horizon flow formalism. The power 
spectrum on large scales is compatible with the WMAP 
7 year data, while the perturbations on small scales are 
sufficiently large that PBHs may be over-produced. The 
Stewart-Lyth calculation is in good agreement with the 
numerical calculation until the final few e-folds of infla- 
tion. On these small scales, the assumptions that are em- 
ployed in the Stewart-Lyth calculation break down, and 
the numerical calculation finds a significant enhancement 
of the amplitude of the perturbations. 



III. RESULTS 

We use the modified flow algorithm described in 
sec. Ill Al to generate a large ensemble (250,000) of infla- 
tion models. In fig. [2] (top row) we plot the cosmological 
observables for all models which are able to sustain the re- 
quired number of e-foldings of inflation, N cos . In around 
7% of the models inflation ends naturally via en = 1 and 
these largely populate the concentrated diagonal feature 
seen in the left hand plots as well as the r = line (c.f. 
Ref. [H|). In the remaining 93% of models, en — > and 
it is assumed, c.f. Ref. Q, that a secondary mechanism, 
such as hybrid inflation, acts to end inflation in these 
cases. Large positive running is in principle allowed (see 
top right plot), however these models may have large 
amplitude perturbations on small scales and hence over- 
produce PBHs. 

To apply the PBH constraints we use the Stewart-Lyth 
expression for the power spectrum, eq. (|15l) . to iden- 
tify inflation models where the amplitude of the pertur- 
bations on small scales which exit the horizon close to 
the end of inflation is large, and may lead to the over- 
production of PBHs. For these models, we then carry 
out an accurate numerical evolution of the primordial 
perturbations, as described in Sec. Ill Bl 

The PBH abundance constraints have recently been 
compiled and updated in Refs. [H, [23|. The resulting 
constraints on the amplitude of the power spectrum are 
typically in the range V-r, < 10 -2 — 10 _1 with some scale 
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dependence [22j . To be conservative we use the constraint 
Vn < 10 _1 . The bottom row of fig. [2] shows the cos- 
mological observables for the models which remain once 
those which over-produce PBHs are excluded. The 7% 
of the original models for which inflation ends naturally 
generally have n s < 1 on all scales and so are unaffected 
by the PBH constraints. Of the remaining models, in 
which inflation continues indefinitely in the absence of a 
secondary mechanism, 92% are excluded by PBH over- 
production. Of the models initially generated, only ap- 
proximately 1% end via a secondary mechanism and do 
not overproduce PBHs. With an accurate numerical cal- 
culation of the perturbations the number of these models 
decreases by approximately 10%. Large positive running 
is now excluded as expected (see bottom-right plot). 

Cosmological constraints on dn s /dlnfc eliminate a sig- 
nificant fraction of the models generated using flow al- 
gorithms [ll|. A full MCMC analysis of cosmological 
data is beyond the scope of this work, however a simple 
application of the observational constraints shows that 
a significant fraction of cosmologically viable models are 
excluded by PBH constraints. Of the models generated 
using our modified flow analysis which have cosmological 
observables within the 3tr ranges found by WMAP7 [12J 
19% are excluded by PBH over-production. This illus- 
trates that in the era of precision cosmological measure- 
ments PBH still provide a powerful constraint on infla- 
tion models. 



IV. CONCLUSIONS 

We have applied constraints on the primordial power 
spectrum from the overproduction of primordial black 
holes to inflation models generated by a modified flow 
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algorithm. The amplitude of inflationary perturbations 
is usually calculated using the Stewart-Lyth [1] expres- 
sion, however for scales which exit the horizon close to 
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of viable models decreases if the power spectrum is cal- 
culated numerically. Of the models generated using our 
modified flow analysis which have cosmological observ- 
ables within the 3tr ranges found by WMAP7 [H 19% 
are excluded by PBH over-production. 

We conclude that PBH constraints provide a signif- 
icant constraint on models of inflation. Furthermore to 
exploit their full power an accurate numerical calculation 
of the amplitude of primordial perturbations on small 
scales, which exit the horizon close to the end of infla- 
tion, is required. 
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